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ABSTRACT1 

This  paper  investigates  a  vibrational  technique  for  the  non  destructive  identification  of 
damage  in  beams.  There  has  been  much  work  investigating  the  variation  of  natural  frequencies 
with  damage,  but  much  less  on  the  effect  of  damage  on  mode  shapes.  Very  little  attention  has 
been  directed  to  off-resonant  performance.  This  paper  develops  a  technique  for  identifying  the 
location  of  structural  damage  in  a  beam  using  experimental  frequency  response  function  data. 
The  procedure  specifically  has  no  requirement  for  knowledge  about  the  undamaged  structure. 
Although  the  technique  is  based  on  modal  superposition,  a  modal  analysis  of  the  measured  data 
is  not  required,  and  offers  no  additional  information. 

The  procedure  is  initially  developed  using  a  finite  element  model.  However,  the  major 
demonstration  in  this  paper  concerns  the  analysis  of  experimental  data  taken  from  a  steel  beam. 
This  was  damaged  with  an  increasingly  deep,  2.8  mm  wide  slot.  It  is  shown  that  the  technique  is 
highly  sensitive,  and  can  locate  minor  damage.  For  example,  the  procedure  correctly  located  the 
slot  when  it  was  only  0.05  mm  deep  (0.8%  reduction  in  thickness).  While  traditionally  the  onset 
of  damage  is  often  indicated  by  a  change  in  natural  frequencies,  at  this  small  level  of  damage  the 
change  in  natural  frequencies  is  within  the  bounds  of  experimental  error,  and  could  not  be  used. 
Other  previously  reported  damage  detection  procedures  were  not  successful  at  locating  the  slot 
until  it  was  more  than  2.25  mm  deep. 


'This  Technical  Report  includes  a  complete  family  of  graphs  and  diagrams,  some  of  which  are  referenced  out  of  order, 
and  some  of  which  are  not  mentioned  at  all  in  the  text.  References  are  in  an  unconventional  format  to  aid  preparation  of 
other  papers,  and  again,  not  all  are  mentioned  in  the  text. 
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1.  INTRODUCTION 

Vibration  measurements  have  been  used  for  many  years  to  identify  damage  in  a  structure.  Some 
of  the  simplest  techniques  monitor  the  vibration  of  running  equipment  (for  example,  aircraft 
engines),  and  if  the  level  exceeds  a  prescribed  quantity  or  condition,  the  monitoring  system  sets 
an  alarm  or  automatically  shuts  down  the  faulty  equipment.  Most  of  these  "engine  health" 
techniques  are  limited  in  that  they  cannot  identify  the  location  of  the  damage,  only  that  damage 
exists.  Increasingly,  it  is  becoming  possible  to  use  vibration  techniques  to  locate  the  damage 
within  a  structure. 

Much  of  the  previously  reported  work  uses  modal  analysis  to  locate  damage  in  trusses  and 
frameworks.  This  work  was  brought  about  by  the  need  to  locate  damage  in  difficult-to-reach 
structures  such  as  the  underwater  part  of  offshore  oil  platforms,  and  space  structures  such  as 
satellites  and  the  space  station.  In  these  structures  the  change  in  natural  frequencies  caused  by 
the  failure  of  a  single  member  can  be  large  enough  to  be  detected  and  used  to  identify  damage. 
Many  damage  detection  techniques,  including  [{Raa-{Raq],  essentially  treat  the  frames  as 
discrete  systems,  and  use  different  methods  of  comparing  modal  behavior  between  the 
undamaged  and  damaged  states. 

Compared  to  trusses  and  space  frames,  there  is  much  less  reported  work  on  beams  and  plates. 

The  main  focus  of  this  work  is  often  on  composite  materials  and  concrete  structures,  such  as 
bridges.  For  example,  Javor  [{Rbb]  gave  international  guidelines  recommending  that  the 
fundamental  frequency  of  bridges  is  monitored  for  long-term  observations.  Miller  et  al.  [{Rbc] 
reported  on  a  reinforced  concrete  bridge,  and  compared  experimental  and  finite  element  results 
to  identify  damage  to  the  shoulders  of  the  bridge  using  differences  in  natural  frequencies,  and 
anomalies  in  mode  shapes.  Liang  et  al  [{Rbw]  presented  an  energy  transfer  method  using 
natural  frequency  and  mode  shape  data  for  detecting  post-earthquake  damage  in  bridges,  and 
Casas  and  Aparicio  [{Rbd]  minimized  a  scalar  performance  error  between  footprint  and 
measured  mode  shapes  and  natural  frequencies. 

The  mode  shapes  of  continuous  structures  change  with  damage,  and  Yuen  [{Rbi]  documented 
the  change  in  the  fundamental  mode  shape  for  a  cantilever  with  respect  to  the  location  of 
damage.  A  displacement  mode  shape  by  itself  is  not  a  reliable  indicator  of  damage  location,  and 
several  authors  have  considered  using  strain  mode  shapes.  Chang  et  al.  [{Rbk]  investigated  the 
sensitivity  of  modal  parameters  to  damage  and  compared  strain  measured  mode  shapes  with 
displacement  mode  shapes.  They  concluded  that  strain  shapes  were  more  effective  at  identifying 
the  location  of  damage  than  displacement  shapes.  This  finding  is  consistent  with  a  study  by 
Pandey  et  al.  [{Rbl],  Curvature  is  proportional  to  the  surface  strain,  and  Pandey  showed  that  the 
difference  between  the  curvature  mode  shapes  for  an  intact  and  damaged  beam  could  locate  a 
localized  change  in  elastic  modulus  of  about  30%. 
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The  relative  sensitivity  of  strain  shapes  to  damage  has  lead  to  a  variety  of  strain  energy  based 
techniques,  including  [{Rbq,  {Rbu],  Damage  is  generally  located  by  comparing  some  strain 
energy  parameter  with  the  same  quantity  for  the  structure  in  an  undamaged  state.  Varying 
degrees  of  success  and  sensitivity  are  reported.  When  an  undamaged  structure  is  not  available 
for  analysis  (or  by  preference  for  some  researchers),  a  validated  finite  element  model  of  the 
undamaged  structure  is  often  used  as  a  "footprint".  However,  requiring  this  footprint  can  limit 
the  application  of  these  techniques.  A  study  by  Ratcliffe  [{Rbq]  introduced  a  technique  whereby 
the  strain  mode  shapes  were  calculated  from  resonant  displacement  mode  shapes.  These  were 
then  processed  to  yield  a  "Difference  Function"  that  located  the  damage.  Data  from  an 
undamaged  structure  were  specifically  not  required. 

One  of  the  limitations  of  all  experimental  modal-based  techniques  is  their  sensitivity  to  data  and 
curve  fit  errors.  Depending  on  the  test  technique  and  the  algorithms  and  curve  fits  used  for  the 
modal  analysis,  there  can  be  an  unexpectedly  large  errors  in  the  modal  constants.  Variations  in 
the  structure  itself  during  testing  can  also  adversely  effect  the  modal  data,  and  surprisingly  there 
is  little  attention  to  these  errors  in  the  literature  [{Rbx].  The  errors  in  modal  constants  are 
particular  significant  since  many  of  the  techniques  rely  on  the  differences  between  the  constants. 

The  technique  presented  here  analyzes  unprocessed  frequency  response  data,  including  the  off- 
resonant  data,  and  thus  eliminates  all  errors  associated  with  a  modal  analysis.  The  procedure  of 
[{Rbv]  is  applied  at  each  frequency  in  the  spectra,  and  a  contour  plot  of  a  Damage  Index  versus 
beam  position  gives  a  strong  visual  indication  of  the  location  of  the  damage.  The  procedure  is 
validated  with  experimental  data  captured  from  a  steel  beam.  A  small  stiffness  reduction  of 
about  2%,  caused  by  locally  reducing  the  thickness  of  a  6  mm  thick  (1/4  inch)  beam  by  0.05  mm 
(2/1000  inch),  gave  virtually  no  measurable  change  in  natural  frequencies.  The  technique 
presented  here  successfully  located  this  small  amount  of  damage.  Other  techniques  presented  in 
the  literature  did  not  locate  the  damage  until  the  thickness  reduction  was  about  ??%. 
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2.  STRAIN  MODE  SHAPES 

Strain  mode  shapes  are  more  sensitive  indicators  of  damage  location  than  displacement  mode 
shapes,  but  compared  to  displacement  shapes,  they  are  often  more  difficult  to  measure  directly. 
Assuming  pure  bending  the  surface  strain  for  a  beam  can  be  related  to  the  radius  of  curvature,  R: 


Especially  for  off-resonant  conditions,  the  displacements  are  small  and  the  radius  of  curvature 
can  be  approximated  by: 
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Experimental  displacement  mode  shape  data  (and  eigenvectors  from  a  finite  element  analysis) 
are  discrete  in  space.  Considering  Equation  (2),  the  strain  mode  shape  can  be  estimated  using  a 
finite  difference  approximation  such  as  Laplace's  Difference  Equation  [(Rbm],  Normally 
Laplace's  equation  is  applied  to  problems  involving  two  dimensions.  In  this  case  the  mode 
shapes  for  a  beam  are  one-dimensional,  and  at  each  resonance  the  Laplacian  of  the  discrete 
mode  shape,  ah  is  calculated  using  Laplace's  Equation  as: 

^  =  Ki  +  a,-i>  '  2ai  (3) 


In  this  formulation  it  is  assumed  that  the  grid  locations  are  equally  spaced  along  the  beam,  and 
that  the  absolute  values  of  the  Laplacian  are  not  of  interest. 
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3.  DAMAGE  DETECTION  WITH  STRAIN  SHAPES 

As  has  been  shown  by  several  authors,  including  [{Rbl  {Rbt  and  {Rbv],  a  strain  mode  shape  by 
itself  can  sometimes  reveal  the  location  of  damage.  As  an  example,  consider  a  uniform  beam 
with  a  localized  reduction  in  stiffness.  This  is  representative  of  crack  damage,  and  may  be 
analyzed  by  reducing  the  thickness  of  one  element  of  the  finite  element  model,  but  leaving  the 
mass  matrix  unchanged  [{Rai]  and  [{Rbl],  Here  the  localized  percentage  reduction  in  bending 
stiffness,  El,  is  called  the  percentage  damage.  The  Laplacian  calculated  using  Equation  (3)  for 
the  first  bending  mode  of  a  finite  element  free-free  beam  with  75%  damage  between  nodes  #7 
and  #8  is  shown  in  Figure  1.  This  level  of  damage  is  severe,  and  is  easily  located.  With  less 
severe  damage  the  Laplacian  retains  its  characteristic  shape,  but  the  effect  is  less  pronounced. 
This  is  shown,  for  15%  damage  to  the  same  beam,  in  Figure  2.  As  described  in  the  Introduction, 
several  techniques  rely  on  the  difference  between  a  footprint  and  damaged  data  in  order  to 
enhance  this  diminishing  irregularity. 


4.  MODIFYING  THE  STRAIN  SHAPES 

As  the  level  of  damage  further  reduces,  the  distinctive  shape  of  the  Laplacian  continues  to 
become  less  pronounced.  There  are  several  methods  that  can  be  used  to  enhance  the  irregularity 
in  the  graph,  such  as  piecewise  linearization  and  cubic  spline.  The  one  progressed  in  this  paper 
was  initially  proposed  in  [{Rbv],  The  method  is  to  fit  a  cubic  polynomial  to  the  Laplacian,  and 
determine  a  Damage  Index  calculated  as  the  numeric  difference  between  the  cubic  and 
Laplacian.  A  separate  cubic  is  determined  for  each  element  of  the  Laplacian  in  turn,  with  the 
coefficients  being  determined  from  the  data  on  either  side  of  the  element,  but  excluding  the 
actual  element.  For  example,  when  the  cubic  calculated  for  the  i*  element  of  the  Laplacian, 
at  position  x;  along  the  beam,  is  defined  as: 

Po  +  PlXi  +  P2Xi  +  P3X?  (4) 


The  coefficients  ph  p2  and  p3  are  determined  using  Laplacian  elements: 


S£a 
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The  Damage  Index,  dh  is  then  calculated  from  the  cubic  and  the  Laplacian. 
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The  Damage  Index  for  15%  damage  between  nodes  #7  and  #8  is  shown  in  Figure  3.  While  the 
'  Damage  Index  does  not  identify  the  location  of  damage  as  accurately  as  the  Laplacian,  it  is  a 
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much  more  sensitive  parameter,  in  that  small  amounts  of  damage  trigger  significant  features  in 
the  graph.  Furthermore,  as  is  shown  later  in  the  experimental  section,  the  procedure  includes  an 
inherent  smoothing  process  that  helps  reduce  the  effects  of  measurement  errors. 


5.  OPERATING  MODE  SHAPES 


Most  previous  experimental  work  has  focused  on  using  resonant  modal  data.  As  discussed 
above,  the  modal  data  are  influenced  by  a  variety  of  errors,  and  since  the  techniques  may  use  the 
differences  between  modal  constants,  the  effectiveness  of  the  damage  detection  procedures  can 
be  seriously  degraded.  The  errors  associated  with  curve  fitting  can  be  eliminated  by  considering 
operating  mode  shapes  determined  directly  from  raw  frequency  response  function  data.  The 
general  response  of  a  structure  subject  to  harmonic  excitation  at  circular  frequency  o>  can  be 
determined  by  a  linear  superposition  of  the  eigenvectors.  For  single-point  harmonic  force 
excitation  at  circular  frequency  o),  at  the  j*  coordinate,  the  spatial  response  V;  at  the  i*  coordinate 
is  given  by: 


f, 


E 


r  =  1 


-  G)2  +  2  / 


(7) 


where  N  is  the  number  of  modes  included  in  the  analysis,  A^y  is  the  modal  constant  associated 
with  the  r*  natural  frequency  and  spatial  coordinates  i  and  j,  and  cor  is  the  r*  natural  frequency. 
Viscous  damping  is  assumed  in  the  equation,  but  is  not  essential  to  this  exposition.  This  well- 
known  equation  shows  that  off-resonant  displacement  shapes  (or  "operating  mode  shapes"  as 
they  may  be  called)  are  a  superposition  of  the  resonant  mode  shapes. 
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6.  DETECTING  DAMAGE  WITH  THE  OPERATING  MODE  SHAPES 


The  Laplacian  can  now  be  determined  for  an  off-resonant  excitation  frequency.  Combining 
equations  (3)  and  (7)  yields  the  general  Laplacian  as  a  function  of  excitation  frequency  and 
discrete  position  i  along  the  beam: 


ST  = 


(v  +  v  ,)  -  2  v 

v  i+1  i- V  i 

Jf  (a  ,  +  a  ,  -  2  a  \ 

V'  \  r,i+ \J  _ MJ/ 

r=  1  |(0*  -  CO2  +  21'C^Wwj 

N  Cg 

_ rj _ 

r  =  1  |(0*  -  0)2  +  2/Cr^rw| 


(8) 


where  Sf;*  indicates  the  non-resonant  Laplacian  function,  and  $£rji  is  the  Laplacian  calculated  at 
resonance  for  the  r*  mode.  This  equation  shows  that  the  Laplacian,  or  curvature  shape,  for  any 
excitation  frequency  is  a  superposition  of  the  Laplacian  functions  calculated  at  resonance.  It  was 
shown  previously  that  damage  causes  a  feature  in  the  resonant  Laplacian  graphs,  and  so  a  feature 
will  also  exist  in  Laplacian  functions  calculated  from  non-resonant  data. 

The  damage  detection  procedure  is  therefore  to  measure  frequency  response  functions  for  an 
equally-spaced  mesh  of  test  coordinates  along  the  length  of  a  beam.  The  data  are  captured  on  a 
digital  spectrum  analyzer,  and  have  a  finite  number  of  frequency  lines.  At  each  frequency  line 
in  the  spectra  the  Laplacian  of  Equation  (3)  is  determined  from  the  operating  mode  shapes 
defined  by  the  unprocessed  frequency  response  data.  The  modulus  of  the  complex  Laplacian  is 
then  graphed  as  a  contour  plot  of  frequency  versus  position  along  the  beam.  In  order  to  improve 
visualization  some  form  of  normalization  is  required.  The  most  effective  procedure  is  to 
normalize  the  Laplacian  at  each  frequency  such  that  the  root  mean  square  value  is  one. 

At  each  frequency  the  complex  Damage  Index  of  Equations  (4)-(6)  is  determined  from  the 
Laplacians,  and  its  modulus  is  normalized  and  graphed  as  a  contour  plot. 

Inspection  of  the  contour  plots  can  reveal  a  number  of  different  conditions.  Remembering  that 
damage  is  located  at  a  specific  position  on  the  beam,  and  that  the  location  of  the  damage  is 
independent  of  frequency,  trends  on  the  contour  plot  that  are  parallel  to  the  frequency  axis  are 
indicative  of  structural  irregularities.  Features  that  are  parallel  to  the  spatial  axis  only  occur  at 
one  frequency,  and  are  associated  with  experimental  error  at  resonance.  Features  that  vary  with 
frequency  show  dependency  on  excitation,  and  are  typically  associated  with  antiresonances. 
Figure  4  shows  the  Laplacian  calculated  from  experimental  frequency  response  data  taken  from 
a  steel  beam  with  no  visually  obvious  damage  (the  full  experimental  setup  and  results  are 
described  below).  The  figure  demonstrates  the  effect  of  both  resonant  and  antiresonant 
behavior. 
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Figure  5  shows  the  Damage  Index  calculated  from  the  Laplacian  data  in  Figure  4.  There  appears 
to  be  a  structural  feature  at  about  coordinate  #27  which  was  not  identified  in  the  Laplacian  of 
Figure  4.  This  feature  is  coincident  with  one  of  the  1/4-inch  suspension  holes.  While  it  is  hoped 
that  the  algorithm  has  detected  the  hole,  the  claim  has  not  been  investigated  fully,  and  is  not 
made  in  this  paper.  Features  on  the  contour  plots  are  best  observed  when  the  plots  are  in  color. 

In  the  reduction  to  monochrome,  the  lowest  valued  contour  was  erased.  All  peaks  above  an 
arbitrary  contour  at  about  RMS=0.35  were  flood-filled  with  black.  Also,  plotting  the  square  of 
the  modulus  generally  produced  an  image  more  amenable  to  monochromatic  rendition,  and  is  the 
procedure  used  for  the  figures  shown  here. 


7.  EXPERIMENTAL  DEMONSTRATION 

7.1  EXPERIMENTAL  CONFIGURATION 

A  flat  steel  beam,  0.9  m  x  0.075  m  x  6.3  mm  thick  (36  inch  x  3  inch  x  0.249  inch),  was 
suspended  through  two  6  mm  holes  positioned  1/4  and  3/4  down  the  length,  and  centered  12  mm 
from  one  edge  of  the  beam.  Suspension  was  with  monofilament  and  rubber  bungee  cords.  A 
uniform  mesh  of  test  coordinates  was  placed  on  the  beam  as  shown  in  Figure  6.  The  beam  was 
initially  tested  by  measuring  frequency  response  functions  between  each  test  coordinate  and  an 
accelerometer  fixed  at  coordinate  #1  (centerline,  left  end).  Modal  data  are  not  required  for  the 
damage  detection  procedure,  however  for  visualization  and  to  assess  data  quality,  the  data  were 
subject  to  a  modal  analysis  using  the  commercial  STAR  software.  For  the  frequency  range 
tested  (0-1.5  kHz),  13  natural  frequencies  were  identified,  there  being  eight  'bending'  and  five 
'torsional'  modes.  Viscous  damping  was  typically  0. 1-0.5%. 

The  beam  was  then  damaged  with  an  increasingly  deep  2.8  mm  (7/64  inch)  wide  groove  milled 
across  the  width  of  one  face  of  the  beam  and  centered  0.33  m  (13  inch)  from  one  end.  The 
width  of  the  groove  was  kept  as  small  as  possible,  commensurate  with  being  able  to  use  a  depth 
gage  to  measure  its  depth.  The  groove  introduced  a  stiffness  change,  with  minimal  effect  on  the 
mass  of  the  beam.  The  beam  was  tested  with  a  total  of  ??  different  groove  depths,  ranging  from 
2/1000  inch  to  ??.  At  each  groove  depth  the  frequency  response  function  data  collection  for  the 
72  coordinate  locations  and  the  modal  analysis  were  repeated.  For  comparison  only  the 
percentage  change  in  experimental  natural  frequencies  for  the  first  four  bending  modes  are 
shown  in  Figure  7. 

7.2  EXPERIMENTAL  RESULTS 

For  all  the  contour  plots,  the  "Beam  Position"  is  the  grid  number.  There  were  36  grids,  starting 
0.5  inch  from  one  end,  with  the  last  grid  being  35.5  inch  from  that  end.  The  groove  is  marked 
on  the  plots  as  a  straight  line  mid  way  between  grids  #12  and  #13.  The  "Frequency  Index"  on 
the  plots  is  the  number  of  the  frequency  line  in  each  spectrum.  Lines  50-700  from  801-line 
0-1.56  kHz  spectra  were  used,  thus  the  actual  frequency  is  calculated  as  (frequency  index 
+49)x  1560/800  Hz. 
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Figure  9  shows  the  Damage  Index  for  a  groove  depth  of  0.002  inch.  Even  at  this  very  low  level 
of  damage  (0.8%  reduction  in  thickness)  the  Damage  Index  shows  a  trend  (an  area  of  black 
parallel  to  the  frequency  axis)  near  beam  position  12-13.  This  is  coincidental  with  the  actual 
location  of  the  damage.  This  feature  in  the  Damage  Index  contour  plot  becomes  further 
enhanced  as  the  groove  is  made  deeper  into  the  beam.  Depths  of  0.005  inch  or  greater  cause  the 
dominant  feature  on  the  Damage  Index  plot  to  correctly  identify  the  damage  location  without 
ambiguity.  Damage  Indices  for  0.005  inch  and  0.030  inch  depths  are  shown  in  Figures  13  and 
23. 

Other  previously  published  damage  location  procedures  discussed  above  considered  modal  data, 
and  inspect  either  the  curvature  or  the  change  in  curvature  caused  by  the  damage.  With  this 
experiment,  curvature  shapes  did  not  give  any  indication  of  damage  until  the  groove  depth  was 
90/1000  inch  (36%  reduction  in  thickness,  74%  localized  stiffness  reduction).  The  curvature 
shape  with  this  condition  for  mode  1  is  shown  in  Figure  38.  The  difference  in  curvature  between 
a  damaged  and  undamaged  structure  never  successfully  located  the  damage.  The  difference  in 
curvature  for  the  maximum  groove  depth  of  1 10/1000  inch  (44%  reduction  in  thickness,  83% 
localized  stiffness  reduction)  for  mode  1  is  shown  in  Figure  45.  Figure  43  shows  the  difference 
in  curvature  plotted  as  a  contour  plot  for  all  frequencies.  Generous  interpretation  may  indicate  a 
feature  close  to  the  groove  (beam  position  12.5),  but  this  is  not  sufficiently  marked  to  warrant  a 
claim  of  damage.  The  reason  for  the  apparent  failure  of  the  curvature  difference  method  was  not 
investigated.  For  comparison,  the  Damage  Index  for  the  110/1000  inch  deep  groove  is  shown  in 
Figure  42. 


8.  DISCUSSION  AND  CONCLUSIONS 

This  paper  has  presented  a  sensitive  technique  that  uses  non-resonant  frequency  response  data  to 
locate  damage  in  an  otherwise  uniform  beam.  The  procedure  consists  of  converting 
displacement  data  to  strain  mode  shapes  by  applying  a  discrete  finite  difference  approximation. 
The  resulting  strain  shapes  are  further  processed  to  reveal  the  location  of  damage.  The 
procedure  is  applied  separately  at  each  frequency  in  the  frequency  response  functions,  and 
presented  as  a  contour  plot  of  (modulus)2.  Features  on  the  contour  plot  which  appear  for  one 
beam  location,  but  at  many  frequencies,  represent  structural  irregularities.  The  procedure 
correctly  located  a  very  small  quantity  of  damage,  2/1000  inch  in  a  250/1000  inch  thick  beam. 

Previous  work  by  several  authors  considered  locating  the  damage  using  the  difference  in 
curvature  of  a  beam  before  and  after  damage.  Curvature  shapes  for  the  fundamental  mode 
located  the  damage  once  the  depth  was  90/1000  inch  (about  one-third  of  the  depth).  Resonant 
Difference  in  curvature  procedures  did  not  locate  the  damage  at  the  maximum  groove  depth  of 
110/1000  inch. 
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Figure  1:  Laplacian  for  50%  damage.  Figure  2:  Laplacian  for  5%  damage. 


Figure  3:  Damage  Index  for  5%  damage 
(free-free). 
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Figure  4:  Laplacian  undamaged  beam 


Figure  5:  Damage  Index  undamaged  beam 
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Figure  8:  Laplacian  2/1000  inch 
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Figure  9:  Damage  Index  2/1000  inch 


Figure  10:  Laplacian  3.5/1000  inch 
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Figure  ll:Damage  Index  3.5/1000  inch 
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Figure  14:  Laplacian  8/1000  inch 
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Figure  15:  Damage  Index  8/1000  inch 
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Figure  16:  Laplacian  14/1000  inch 
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Figure  17:  Damage  Index  14/1000  inch 
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Figure  18:  Laplacian  19/1000  inch 
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Figure  19:  Damage  Index  19/1000  inch 
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Figure  20:  Laplacian  24/1000  inch 
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Figure  21:  Damage  Index  24/1000  inch 
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Figure  22:  Laplacian  30/1000  inch 
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Figure  23:  Damage  Index  30/1000  inch 


VIWJ 

IW1U 


m 


[■■■S 

rnilWBIH 

giniiinniii 
giiiiiiiiiinvil 
EiiiiiiniiswSI 


SimTllillflTIllllV 


BEAM  POSITION 

BEAM  POSITION 

Figure  24:  Difference  in  Curvature  for 
30/1000  inch.  Mode  1. 

Figure  25:  Difference  in  Curvature  for 
30/1000  inch,  Mode  2. 
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Figure  26:  Laplacian  difference  30/1000 
inch 
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Figure  27:  Laplacian  39/1000  inch 
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Figure  29:  Laplacian  50/1000  inch 


5  0  15  20  25  30  35 


BEAM  POSITION 

Figure  28:  Damage  Index  39/1000  inch 
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Figure  30:  Damage  Index  50/1000  inch 
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Figure  35:  Laplacian  90/1000  inch 


Figure  36:  Damage  Index  90/1000  inch 


Figure  37:  Laplacian  difference  90/1000 
inch 
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Figure  38:  Laplacian  for  mode  1,  90/1000 
inch 


Figure  39:  Difference  in  curvature  for  mode 
1,  90/1000  inch 
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Figure  40:  Laplacian  for  mode  2,  90/1000 
inch 
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Figure  41:  Laplacian  110/1000  inch 
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Figure  42:  Damage  Index  110/1000  inch 
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Figure  43  :  Laplacian  difference  110/1000 
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Figure  45:  Curvature  difference  mode  1  Figure  44:  Laplacian  Mode  1  110/1000  inch 

110/1000  inch 


Figure  46:  Laplacian  Mode  2  1 10/1000  inch 


Figure  47:  Curvature  difference  mode  2 
110/1000  inch 
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